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£ . Abstract 
i— » 

OO '. We calculate 7/(0), the damping rate for longitudinal on-shell gluons with 

. zero momentum in hot QCD using the hard-thermal-loop (htl) scheme. We 

■ find it to be divergent in the infrared, which means that in this scheme 7/(0) 
is different from 7t(0), the corresponding damping rate for transverse gluons 

(--^ ' which is known to be finite. This result suggests that the htl scheme is infrared 



sensitive and thus may need to be improved upon in this sector. We discuss 



this issue after we present our calculation. 
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Besides their importance regarding the stability of the quark-gluon plasma, gluon damp- 
ing rates have been crucial in better understanding QCD at high temperature T. In early 
works, it has been noticed that in this regime, the determination of the dispersion laws 
for quarks and gluons beyond lowest order using standard perturbation theory is gauge- 
dependent This problem has been emphasized in further works in which the gluon 
damping rates have been calculated to one-loop order in various gauges and schemes and 
different results have been obtained ||. It was then realized that the problem was related 
to the way the expansion in powers of g, the perturbative QCD coupling constant, was per- 
formed: at high T (the hard scale), when the external momenta are soft, i.e., of magnitude 
gT, the standard loop expansion is not anymore an expansion in powers of g @ . It was sub- 
sequently developed an effective perturbative expansion in the framework of a resummation 
scheme of the so-called hard thermal loops (htl) Q]. Using this scheme, the transverse-gluon 
damping rate 7t(0) with zero momentum was shown to be Coulomb-and-covariant-gauge- 
invariant and determined in the strict Coulomb gauge to be finite and positive ||. Later, a 
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generating-functional formalism in the htl approximation was developed || and a relation 
to the eikonal of a Chern-Simons gauge theory was made 0. From there a hydrodynamic 
approach showed that the htl approximation was essentially 'classical' 

However, the htl resummation scheme may be not 'complete' yet in describing the whole 
picture of QCD at high temperature. Indeed, it discusses only the two scales T (hard) and 
gT (soft) whereas with g and T, one has (at least) a hierarchy of scales g n T, n an integer, 
positive if we believe that T is the highest scale in the problem. It may be unnecessary to 
reorganize perturbation theory taking into account this general structure of scales, but there 
are indications that at least the next smaller one g 2 T may be important. Indeed, it is true 
that in the htl scheme, both on-shell longitudinal and transverse gluons acquire a thermal 
mass m g of order gT, but static magnetic fields are not screened yet at this scale [f|; they 
are expected to get so at the hence-called 'magnetic scale' g 2 T. Furthermore, the gluon 
damping rates ^{p) are to lowest order of magnitude g 2 T and it has been noticed that these 
damping rates, when calculated in the htl-resummation scheme at soft but nonzero momenta 
p, exhibit a Ing behavior accompanied with coefficients (T(J. This logarithmic behavior is 
not present in the expression of 7t(0), the limit p — > of 7j(p), and no explanation has 
been provided yet for this discrepancy. These remarks tend to indicate therefore that the 
htl scheme may be sensitive in the 'infrared'. 

We think a more direct indication of the infrared sensitivity of the htl-resummation 
scheme is the observation we made in [JTTJ] that, when calculated solely in this scheme, 
the longitudinal-gluon damping rate at zero momentum 7;(0) may be potentially infrared- 
divergent. Indeed, we have determined in that work the analytic expression to leading order 
of 7/(0). We have obtained for it the following form: 

fn\ 9 2 N C T \ (i) (2)1 , s 

7/(0) = -^r [ a (o + a io\ ; I 1 ) 

N c is the number of colors, dj^ is given in eq (|3|) below; it is a finite number found in || to 
be equal to 6.63538 and is such that the damping rate 7t(0) for transverse gluons with zero 
momentum is equal to 9 „. c |J. The other contribution dj^ is given in eq (^) and 



it is this part of 7/(0) we remarked it contains terms which are divergent in the infrared [3 I 



From a physical standpoint and independently of any calculation scheme, there is no 
difference at zero momentum between longitudinal and transverse gluons, and hence one 
expects any consistent scheme to yield the same result for both damping rates ||. In virtue 
of this physical argument, once 7t(0) was found in the htl scheme to be finite and positive, 
7/(0) was thought to be so and to the best of our knowledge, no explicit calculation of the 
latter in the htl scheme has been reported prior to the one in [11] . With this in mind and 



expecting the htl-resummation scheme to be 'complete', we argued in fly] that the terms 

i {2) 
Ho j 



(2) 

contributing to a l0 , in particular the infrared-divergent ones, when put together, may cancel 
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one another. We show in this letter that this is not the case: when all the expressions are 

(2) 

treated with care, not only we find the part a lQ nonvanishing but the infrared-divergent 
piece survives too, see eq In the sequel, we present first our calculation and we discuss 
the result afterwards. 

The longitudinal gluon damping rate in the strict Coulomb gauge is obtained to lowest 
order by the relation: 



li(p) 



Im *!!;(— iuo,p) 



LL>=u>i{p)+iO+ 



(2) 



where 511; (*H) is the longitudinal-gluon htl (effective) self-energy. u>i(p) is its on-shell 
lowest-order energy given for soft momenta in units of m g by eq (|8|) below. We follow closely 
the notation of fill to which this work is a follow-up. Since the denominator in (|D is given 
in units of m g by d/doj 8H(— iu, p)| w =wi(p) = —2p 2 + ... , to get ji(p = 0), we have to expand 
the imaginary part of the effective self-energy to order p 2 . This calculation is reported in 
TT| and we obtain the result ([!]) where the expressions of a\l and d$ are the following: 
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and: 
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The notation is as follows: @i = Q(k 2 — uj\) where G is the step function and 8\ = 8(k 2 
d k = <9/<9fc etc and pu is a short notation for pi(ui,k), i — 1,2 and the same for p t j. The 
spectral densities p^/ are given by [pUTS]: 



Pt> i(u, k) = h)l (k) [8 (uo - u t;l (k)) -8(lu + u t;l (k))} + Pu(uj, k) B(k 2 - to 2 
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an expression in which the residues }t,i(k) are given by: 



and the cut functions @t,i(u, k) by: 
Pt(u,k) 
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ujt,i{p) are the on-shell energies of the gluon to lowest order. For soft gluons and in units of 



m g , we have: 
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As we said, finite number found in || to be equal to 6.63538. We therefore have 

(2) 

to calculate a l0 . The sort of difficulties one encounters when dealing with the expressions 
involved in eq (f|) can partly be displayed in the following simple example. Consider the 
integral: 
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the residue piece of a term that intervenes in the expression of a l0 . Here e{u) is the sign 
function. We may choose to write 8{k 2 — uj 2 ) = ^rS(k — \u>\) , in which case we obtain: 
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We have used the fact that u>i(k) > 1 for all k and an expansion of $i(h) and ui(k) for small 
k. In the above equation, r\ is an infrared cutoff; 77 <C 1. Thus I is divergent in the infrared. 



But we may also choose to write 5{k 2 — u 2 ) = ^S(k — \uj\) , in which case we have: 
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a divergent result too, but different from that of eq (pUp. Which result should then one 
consider? The ambiguity comes from the fact that the condition k + 1 — uii(k) = is 
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satisfied only at k = 0; it cannot be 'approached from above', so to speak. Thus if one 
chooses to regularize I by defining: 

-oo r+oo j, , r+°° 



l> + QO n+OO J p + OO 7 
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with 1 but kept different from zero until all intermediary steps are performed, for both 
choices of eqs fllUD and ([□]) one obtains: 

I, = 0, (13) 

and the ambiguity is lifted. This is the regularization procedure we adopt in this work. The 
calculation becomes then mainly a matter of straightforwardly disentangling the infrared- 
divergent pieces from the finite ones. It very often necessitates an expansion for small k of 
the residue and cut functions using their definitions given in (j^j and ([?[) respectively and 
some of their derivatives. Also, care must be taken when handling the delta functions that 
occur, especially with their first and second derivatives. Finally, extra care must be given to 
the order in which the integrals over k and u are performed; it can be quite subtle in some 
cases. 

(2) 

We now give the results for the different contributions to a l0 . To ease the notation, we 
denote in a compact way: 

V = 9 f °° dk f°° f°° *± 5(1 - a* - u 2 ) . (14) 

J Tj J — OO 

We first have: 

/O Q /*+00 -i r) 7 

2a*) 0! = - JL _ 3 J dk Y^]7 A(1 " K fc) • (15) 

Next we have: 

/r / k 2 \ 2k 
Vd k [k 3 (l + -\ p lxPn + — {k 2 {l - 3k 2 ) + (1 - Ak 2 )^^ - u^l) PtlPt2 

It is worth mentioning at this stage that all contributions from k — > +oo vanish, indeed as 
it should be since only soft contributions are to contribute. The next piece contains a term 
the residue part of which we have discussed in the example above. It reads: 
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The next pieces are more delicate to handle: they need to work out derivatives of delta 
functions in the presence of other delta functions with different arguments. Nevertheless 
that can be done and we get: 
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The next term is worked out along the same lines as the previous one. We have: 
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The next term is the most tedious: it necessitates the additional expansion of d k /3 t {cu, k) 
with uj = 1 — uj t {k) to order k 2 . When this is done and the contributions put together, we 
get: 
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The last term isn't more difficult. It reads: 
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Putting all these contributions together we obtain: 
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Note that the integral in the above equation (and all the other similar ones actually) is 
finite. We can evaluate it numerically and we get: 



«# = °**» -4.45366. 

tu ^2 



(23) 



Hence, as announced, a$ is nonzero and quite divergent. Via eq (|l|), this means ji (0) is 
itself infrared-divergent. Indeed, using the value a[V = 6.63538, we get: 
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We resume our initial discussion. Remember that the calculation of both j t (0) and 7/ (0) 
is performed solely within the framework of the htl-resummation scheme. This latter treats 
only the two scales T and gT and we have already remarked that the scale g 2 T may play an 
important role in QCD at high T. This suggests that the scheme as it is may not be robust 
enough at scales lower than gT. One natural way this lack of robustness may manifest itself 
is via an infrared sensitivity of the scheme when it comes to deal with quantities that to 
lowest order are already of magnitude g n T, n > 2. We think the divergence in 7/ (0) we 
report in this letter is a direct manifestation of this infrared sensitivity. 

One may then ask why the result is finite for ^ t (0) whereas we find an infrared divergence 
for '-fi (0). It is pertinent to note in this respect that there is a main difference between the 
calculation in the htl scheme of 7/ (0) and that of j t (0): in order to get 7/ (0), we already 
need an unavoidable expansion to order p 2 of the imaginary part of the effective self-energy 
whereas for 7j(0) we do not p|,[TT||. As a matter of fact, we have determined the analytic 
expression of the transverse-gluon damping rate 7t(p) to order p 2 . The coefficient ato, that of 
the zeroth order in p 2 , is finite and as we said, equal to of eq (^) p|. But a preliminary 
investigation shows that the coefficient an, that of order p 2 , has contributing terms which 



are divergent in the infrared . It seems then that the expansion of the gluon effective self- 
energies *I1 of the htl scheme in powers of the soft momentum p is infrared sensitive beyond 
zeroth order. 7/(0) necessitating such an expansion to order p 2 , hence an infrared-divergent 
result. 

One may jump to conclude that the htl scheme is robust in the infrared, and that all 
this simply means that the expansion in powers of the soft momentum of the gluon effective 
self-energies is not valid in the first place. But remember that this expansion is necessary 
for 7/(0); it is 'part of its definition, so to speak. Since on physical grounds the damping 
rates have to be finite for soft momenta, any consistent calculation scheme must allow for 
a small-p expansion with finite coefficients. If the expansion yields infinite coefficients, it 
only means that the scheme in which the calculation is performed needs to be cured and/or 
improved upon, not that the principle of the expansion itself is invalid. It is our opinion 
that this is the case for the htl scheme in the infrared sector. In other words, an infrared- 
improved resummation scheme would still allow for the same expansion but would remove 
the infrared sensitivity of the actual htl scheme. 

In fact, the infrared sector is not the first instance in which the htl scheme yields infini- 
ties. Indeed, colinear singularities do appear when external light-like momenta are involved 
and in response to this, the scheme has been subjected to improvement in [13|] where an 
improved action that removes these singularities is proposed. We think the same treatment 
is necessary in the infrared sector but the situation is more problematic here. Indeed, such 
an improvement would very likely necessitate the determination of the magnetic mass which 
is mostly believed to intervene in a nonperturbative way. It may also be that the story does 
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not end at the scale g 2 T and new physics may manifest itself at lower scales g n T with n > 2, 
though nothing is suggesting this for the moment. In any case, it is pertinent to note that 
besides the result we obtain for 7z(0), what we think is also interesting in this work is that it 
sets a calculable framework in which the infrared divergences appear explicitly and by the 
same token where to test an improved htl scheme when found. In the meantime, we can 
have a very rough estimate of the scale at which the htl scheme starts to be sensitive. If we 
demand that 7;(0) be equal to 7t(0), we may set a^' equal to zero, which yields rj = 0.145 
(in units of m g ). However, one should be prudent in considering this value as an estimation, 
even rough, of the magnetic mass. 

Another interesting point to discuss is the logarithmic behavior in j(p) with p soft but 



nonzero we mentioned in our introductory remarks pOfl . More precisely, one argues that 



lip) ~ fip) l n ( m ma g / m eiec), where f(p) is a well behaved function and the ratio of the 
magnetic mass to the electric one is (at least) of order g. Now as we suggested, our cut-off 
rj can be thought of representing the magnetic scale, and thus in units of m e \ ec (denoted in 
our work by m g ), it would also be of order g at least. We notice then the absence of In rj in 
the divergent behavior we obtain for 7/(0), see ( p4]) . This fact can be used to put back into 
question the validity of the expansion we make. In particular, one may argue that a term of 
the form f(p) Inp should be present in the expansion of the effective self-energies *I1, a term 
more in line with the logarithmic behavior. First of all, such a sort of terms is ruled out 
by the nature of the small— p calculation itself and there is really no subtlety in this as one 
can straightforwardly get convinced by following the steps reported in Also, we stress 
that the htl scheme must allow for the expansion in powers of p and it actually does. What 
simply happens is that the coefficients of the expansion it yields are infrared sensitive, which 
means that the scheme is still yet not 'complete' in this infrared sector; in other words, it 
does not resum diagrams of order g 2 T and smaller whereas our result suggests there is a 
need for such a further resummation in one way or another. Once the scheme is properly 
improved upon, then the new coefficients will cease to diverge. 

But the point about the absence of In t] in the expression of 7;(0) still remains. Indeed, 
if the In rj is present in the expression of 7(p), then one should expect it to remain present 
in that of 7(0). It is not difficult to understand its absence in the transverse case. Indeed, 
when we expand 7*(p) to order p 2 , the coefficient a t o is finite but the coefficient a t i is 
infrared- divergent, in particular the presence of In r\ in it is not ruled out from the outset. 
Therefore, when p is small but different from zero, the coefficients containing the infrared 
pieces dominate over the finite ones like a t Q. But precisely at p = 0, only the coefficient 
a t o survives and it is finite, hence the absence of In rj or any other infrared- sensitive piece. 
However, the situation is more delicate in the longitudinal case: this argument cannot be 
carried forward there since already aio is divergent and does not contain a In rj term. Actually 
the longitudinal situation puts into light a problem we have not raised thus far. Indeed, 
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Besides In 77, we expect to find power-like divergences (some higher than 1/rf) in an- These 
divergences are more severe than In 77 and yet have not been reported in the literature. This 
issue is well beyond the scope of the present work and we defer its discussion until we finish 
the disentanglement of all the infrared pieces in 7t(p) to order p 2 fll4|. 

Finally, we have limited ourselves thus far to discussing the gluonic sector only but 
certainly quarks play an important role too in the structure of hot QCD. They themselves 
acquire a thermal mass mj to order gT and their damping rates start also at the scale g 2 T, 
which means they are not a priori shunned from infrared sensitivity. 
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